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Abstract — The conversion of geometric as well as ellipsoidal heights from GNSS observations to practical heights 
for engineering constructions has necessitated the determination of the local geoid model of areas. Benin City is a 
developing area which requires a local geoid model for conversion of geometric heights to orthometric heights for 
physical developments in the area. This paper is on the best local geoid model of Benin City, Nigeria by comparing 
three gravimetric-geometric geoid models of the study area. GNSS and gravimetric observations were carried out on 
49 points to respectively obtain their coordinates and absolute gravity values. The theoretical gravity values of the 
points were computed on the Clarke 1880 ellipsoid, subtracted from the absolute gravity values and corrected for the 
air (free air) to obtain the free air gravity anomalies of the points. The computed free air gravity anomalies were 
applied to compute the geoid heights of the points using the integration of the modified Stokes integral. Three 
geometric geoid surfaces (plane, second degree and third degree surfaces) were fitted to the computed gravimetric 
geoid heights using the least squares technique to obtain the gravimetric-geometric geoid models of the study area. 
The RMSE of the three gravimetric-geometric geoid models were computed to determine their (the models) accuracy. 
The three gravimetric-geometric geoid models were compared using their accuracy to obtain the most suitable geoid 
model of the study area. The results of the comparison showed that the third degree gravimetric-geometric geoid 
model is most suitable for application in the study area. It is recommended that ellipsoidal heights obtained from 
GNSS observation within Benin City, Nigeria should be converted to orthometric heights using the third degree geoid 
model. 
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I. INTRODUCTION 

The Geoid is one of the Earth’s shapes. Geoid surface 
is used to approximate the physical shape of the Earth. It is 
the equipotential surface of the Earth’s gravity field which 
more or less coincides with mean sea level (Borge, 2013). 
Ubajekwe (2011) also defined the geoid as the 
equipotential surface of the earth’s attraction and rotation 
which coincides on average with the mean sea level in the 
open Ocean. It is the surface which coincides with the 
mean sea level assuming that the sea was free to flow 
under the land in small frictionless channels. Ubajekwe 
(2011) explained further that the mean sea level is not 
quite an equipotential surface owing to non-gravitational 
forces (such as ocean currents, winds and barometric 
pressure variation). In geodetic surveying, the computation 
of the geodetic coordinates of points is commonly 
performed on a reference ellipsoid closely approximating 
the size and shape of the earth in the area of the survey. 
The actual measurements made on the surface of the earth 
with certain instruments are however referred to the geoid. 


Moritz and Hofmann (2005) stated that the geoid coincides 
with that surface to which the oceans would conform over 
the entire earth if free to adjust to the combined effect of 
the earth mass attraction (gravitation) and the centrifugal 
force of the earth’s rotation. They also explained that the 
geoid is a surface along which the gravity potential is 
everywhere equal and to which the direction of gravity is 
always perpendicular when optical instruments containing 
gravity reference levelling devices are properly adjusted 
during observation coincides with the direction of gravity 
and are therefore perpendicular to the geoid. The angle 
between the plumb line which is perpendicular to the geoid 
(the vertical) and the perpendicular to the ellipsoid (the 
ellipsoidal normal) is defined as the deflection of the 
vertical. Civil engineers use the geoid as the reference 
surface for elevations while oceanographers use it for 
studies of ocean circulation, currents and tides. It is also 
valuable to geophysicists for displacement studies, 
geophysical interpretation of the Earth's crust, and 
prospecting (Borge, 2013). 
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The GNSS ellipsoidal height is a geometric height 
which is obtained with reference to the ellipsoid. The 
ellipsoidal height is not suitable for engineering work as it 
is not a practical height determined with respect to the 
geoid. Nowadays, most geodetic observations are carried 
out with GNSS receivers. The spirit levelling, on the other 
hand, is time-consuming and applying it over long 
distances may reduce the reliability of the measurements. 
But there is a relation, as well as a link between the 
ellipsoidal height and the orthometric height. The link 
between the ellipsoidal and the orthometric heights is the 
geoid height. In other words, if the geoid and the 
ellipsoidal heights of a point are known, the orthometric 
height of the point is computed using the relation given by 
Eteje et al. (2018) and Oluyori et al. (2018) as: 

h = H + N (1) 

Where, 

h = Ellipsoidal height; H = Orthometric height; N = 

Geoid height 

The relationship between the heights is further 
explained using the three reference surfaces (the ellipsoid, 
the geoid and the earth surface). The ellipsoidal height is 
determined with respect to the reference ellipsoid; the 
orthometric height is obtained with respect to the geoid 
while the observations are carried out on the earth surface. 
The ellipsoidal height is the vertical distance measured 
from the surface of the ellipsoid along the normal to the 
observation point on the earth surface while the 
orthometric height is the vertical distance measured from 
the geoid along the vertical, as well as the plumbline 
direction to the observation point on the earth surface. 
Figure 1 shows the relationship between the three 
reference surfaces, as well as the ellipsoidal, orthometric 
and the geoid heights (Eteje et al, 2018) 



Fig. 1: Relationship between Orthometric, Geoid and 
Ellipsoidal Heights 
Source: Eteje et al (2018) 

The geoid can be determined using various methods such 
as the gravimetric, geometric, Astro-geodetic, 
transformation and gravimetric-geometric methods 
amongst others. 


The gravimetric method can be carried out by the well- 
known Stokes-integral, equation (2) and the use of 
accurately determined absolute gravity data (Heiskanen 
and Moritz, 1967, Eteje, 2015 and Eteje et al, 2018). 

N = jj SgS(y/)dcr ( 2 ) 

• (7 

Where N is geoid undulation, A g is gravity anomaly, 

s(r) is stokes function, y is normal gravity on the 
reference ellipsoid and R is mean radius of the earth. The 
geometric method is to use the known “geoid heights” at 
some points, which are derived from co-located GNSS- 
determined heights and levelled heights to interpolate the 
geoid heights at other points (Chen, and Luo, 2004). The 
interpolation of the geoid heights at any other point 
involves the use of interpolation models such as bicubic 
model and other models like multiquadratic model, etc. In 
Astro-geodetic method of geoid determination, the geoid 
heights of points are determined with reference to a 
geodetic (reference) station whose geoid height is known. 
The geoid heights differences between points are 
determined using the components of deflection of the 
gravity vector which can be obtained by carrying out 
astronomical and geodetic observations. The astronomical 
observation is carried out to determine the astronomical 
coordinates (astronomical latitude, O and astronomical 
longitude, A ) by observing stars. The geodetic observation 
is used to determine the geodetic latitude, cp , geodetic 

longitude, X and ellipsoidal heights, h as well as the 
azimuths, (X and geodetic distances, i between network 
points. Using the astronomical and geodetic coordinates, 
the components of deflection of the gravity vector can be 
computed. The transformation method involves the use of 
the well-known Euclidean similarity transformation model 
which is used to convert Cartesian coordinates between 
two geodetic reference frames that generally differ in 
terms of three translation parameters (tx, ty, tz), three 
orientation parameters (sx, sy, sz) and a factor of uniform 
spatial scale change (5s) (Eteje et al , 2018). 

The gravimetric-geometric method has to do with the 
use of gravimetrically and geometrically obtained data. 
That is the combination of the gravimetric and the 
geometric methods. Here, the geoid heights of selected 
points are determined with the gravimetric method and a 
geometric surface is fitted to the gravimetric geoid heights 
to enable the interpolation of geoid heights of points of 
interest. Usually, two or more geometric surfaces are fitted 
to the geoid heights of the points and the Root Mean 
Square Error (RMSE), as well as the standard error of each 
of the models, is computed and compared to obtain the 
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most suitable model for application in the study area. The 
fitting of the geometric surfaces to the gravimetric geoid 
heights requires the computation of the model parameters, 
as well as coefficients using the least squares technique. 
The computed model parameters and the coordinates of the 
points are used to develop a program for the 
implementation of the model in the study area. The RMSE 
and as well as the reliabilities of the models are obtained 
by comparing the results from physical measurements 
carried out on some selected points with their 
corresponding results from the determined geoid models. 

Benin City, Nigeria is the capital of Edo State. It is a 
developing area where massive roads constructions are 
presently ongoing. As stated earlier, most observations are 
being carried out using GNSS receivers. Some of the 
preliminary surveys used for the design of the roads were 
carried out using Real-Time Kinematic GNSS observation 
technique. The heights obtained from the observations 
were ellipsoidal heights and they were used without 
converting them to orthometric heights. This is because 
there has not been a local geoid model of Benin City, 
Nigeria. Orthometric heights are required in engineering 
constructions to decide on suitable gradients to direct the 
flow of water at various levels of proposed constructions. 


With orthometric height water can flow from one point to 
another while with ellipsoidal height, it is not possible. 
Today, the Benin City, Nigeria is facing a serious flood 
problem. Several attempts have been made to solve this 
problem but there has not been headway as the height 
differences between points relative to the geoid as the 
mean sea level has not been properly established. This 
study determines and compares three gravimetric- 
geometric geoid models to decide on the best local geoid 
model of Benin City, Nigeria. 

1.1 The Study Area 

Benin City, Nigeria is a City approximately 40 
kilometres north of the Benin River. The City is also 
linked by roads to Asaba, Sapele, Siluko, Okene, and 
Ubiaja and is served by air and the Niger River delta ports 
of Koko and Sapele. The City is made up of three Local 
Government Areas, Oredo LG A, Ikpoba Okha LG A and 
Egor LGA. It has a total population of 1,749,316 according 
to 2019 NPC projection. It covers a total area of about 
1,204 km 2 . Benin City is bounded by UTM zone 31 
coordinates 660000mN and 712500mN, and 770000mE 
and 815000mE. Ligures 2a and b show the maps of the 
study area. 



Fig. 2 a: Map of Edo State 


Fig. 2b: Map of Benin City 
Source: Ministry of Lands and Surveys, Benin City 


1.2 Normal Gravity Computation 

The normal, as well as the latitude gravity, is computed on a specified ellipsoid. Eteje et al. (2018) gave the model for the 
computation of the theoretical gravity on 
the Clarke 1880 ellipsoid as 
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Where, 


8rciarke\SS0 ( B ) 


9.78051938 


(l + 0.0052474<36sin 2 cp' 
(-0.0000087SB5sin 2 2 rp^ 


ms 


8t ciarke\ wo = Theoretical gravity on the Clarke 1880 ellipsoid 

cp = Station latitude 


(3) 


1.3 Gravity Anomaly 

The gravity anomaly, Ag, is the difference between 
the observed gravity value (g) reduced to the geoid, and a 
normal, or theoretical, computed gravity value (/ Q ) at 

the mean earth ellipsoid, where, the actual gravity 
potential on the geoid equal the normal gravity potential 
at the ellipsoid, at the projection of the same terrain point 
on the geoid and the ellipsoid respectively, that is 
(Dawod, 1998 and Eteje et al ., 2019) 


to Aziz et al (2010), this correction is added to the 
observed gravity because the increased radial distance of 
the station from the centre of the Earth results in a lower 
observed gravity value than if the station were at the local 
datum. The formula to calculate the magnitude of the 
reduction in practice is given by Eteje et al. (2019) as 

8 fa =~ — H s = -308.6HpGal 

r (5) 

= -0.3086 HmGal 


A 8 = 8 ~Yo ( 4 ) 

Considering the nature of the topography of the earth 
surface, which is irregular in shape, there are two basic 
types of gravity anomalies (free air and Bouguer 
anomalies). In this study, it was only the free air 
correction that was applied. That is free air gravity 
anomalies that were used. 

1.4 Free Air Correction 

This is the first step for reducing topography effects. It 
simply corrects for the change in the elevation of the 
gravity meter, considering only air (hence a free-air) 
being between the meter and selected datum. According 


N = 


rAg 


8 r 


-6 sin 2 y/ o ln<! sin 


f W_P 

v2y 


■ 24sin ; 


Wo 


12sin 


+ sin 


Wo 


Where, 

H = Station orthometric height in metres 
g = Mean value of gravity (980500 mGal) 
r = Mean radius of the Earth 

1.5 Integration of Stokes’s Formula 

According to Eteje et al (2018), using the modified 
Stokes integral given in equation (2), the geoid heights of 
points are computed if their gravity anomalies and 
geographic coordinates are known. Featherstone and 
Olliver (1997) gave the integration of equation (2) as 


/ \ 

l 

/ \ 


/ \ 

\ 

Wo 

> + 16sin 

Wo 

+ 12sin 2 

Wo 


2 J 

J 

U J 


l 2 J 



-Acosy/ o +5cos2^ o -1 

) 


Where N is the geoidal height of individual point, y/ 0 is 
the surface spherical radius, y is the theoretical as well 
as normal gravity, A g is the gravity anomaly and r = R 
is the mean radius of the earth. 

1.6 Computation of Surface Spherical Radius, l// () 

The surface spherical radius, y/ 0 is computed as 
(Shrivastava et al, 2015) 

cosy/ = sin^sin^ 1 + cos^cos^ 1 cos(A l - A) (7) 
Where, 

cp — Mean latitude of the points 
cp 1 = Latitude of individual point 
A — Mean longitude of the points 


A 1 = Longitude of individual point 

1.7 Computation of Combined Topographic Effect 

To obtain a precise geoid height of a point, the 
combined topographic effect is calculated and applied to 
the computed geoid height of the point. The formula for 
the computation of the combined topographic effect, 

SN T c Z b is given as (Sjoberg, 2000 and Kuczynska- 


Siehien et al., 2016): 


SN Topo 

Uly Comb 


IjtGp 


Y L 


H 2 + — H 2 
3 R 


(8) 


where G is the earth gravitational constant, p is density, 

R is the mean radius of the earth and H is the orthometric 
height of observation point which can be obtained from 
the DTM of the area. 
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1.8 Geometric Geoid Surfaces 

According to Eteje et al. (2018), geometric geoid 
surfaces are mathematical interpolation surfaces fitted to 
geoid heights to enable geoid heights of new points to be 
determined using variable such as geographic or 
rectangular coordinates of the points. These surfaces 
include plane surface, bi-linear surface, second degree 
surface, third degree polynomial and fifth degree 
polynomial. Eteje et al. (2018) further explained that the 
surface to be adopted as well as the degree and order of 
the polynomial depends on the size of the study area and 
the variation of the geoid heights. For a small area, the 
plane surface is used, for a relatively large area, the 
second and third order polynomial surfaces are used. 

1.9 Plane Surface 

The plane surface as given by Odumosu et al (2016) is 
N = (h- H) = N 0 + N x e + N 2 n + S N (9) 


The least squares solution to equation (9) can be further 
mathematically expressed as (Odumosu et al , 2016): 


(h- 

H) = 

-N a 

+ N 

{ e + N 2 n 

(10) 

C H ) 

= h - 

~N 0 

-N. 

{ e-N 2 n 

(11) 





( N ' 1 

iy 0 



(h~ 

H) = 

= d 

e ri) 

N 1 


(12) 





K N 2 y 




Where X is the required parameters = 


N 0 

N\ 

KN 2 ; 


, e and n are 


respectively easting and northing coordinates of the 
selected points. 

1.10 Polynomial Surface 

The polynomial surface for local geoid heights 
interpolation given by Kirici and Sisman (2017) and 
Oluyori et al. (2018) is 


N (*,y)=1L 


(13) 


1=0 j=k-i 
i =0 


Where, 

a tj = Polynomial coefficients 

m = Degree of polynomial 
x, y = Plan coordinates of points 

In applying the polynomial, the degree should be chosen 
and the polynomial should be formed for the chosen 
degree. Kirici and Sisman (2017) gave the third degree 
polynomial as 


N — (2qq + ci^qX + cLq^Y + U2qX + ct^XY + ct^^X -\-ci 21 X Y + ct^XY -\-cIq?Y^ (14) 

Also, Sanlioglu et al (2009) gave the second degree polynomial, Multi-quadratic model as: 

N = a 0 + a x X + a 2 Y + a 3 X -\-ct^Y + ct^XY + ct^X Y-\-ctjXY -\-ct%X Y 15) 


Where, 

Y = ABS(y-y 0 ) 

X = ABS (x-x o ) 

y = Northing coordinate of observed station 
X = Easting coordinate of observed station 
y o = Northing coordinate of the origin (average of 
the northing coordinates) 

X 0 = Easting coordinate of the origin (average of the 
easting coordinates) 

1.11 Observation Equation Method of Least Squares 
Adjustment 

The fitting of the polynomial interpolation surface to a 
set of geoid heights requires the model parameters 
(variable coefficients) to be computed. The computation of 
these coefficients is done by observation equation method 
of least squares adjustment technique. The functional 
relationship between adjusted observations and the 


adjusted parameters as given by Eteje and Oduyebo (2018) 
is: 

L a =F(X a ) (16) 

Where L a = adjusted observations and X a = adjusted 

parameters. Equation (16) is a linear function and the 
general observation equation model was obtained. The 
system of observation equations is presented by matrix 
notation as (Ono et al., 2018): 

V = AX-L (17) 

Where, 

A = Design Matrix, 

X = Vector of Unknowns 
L = Observation Matrix. 

V = Residual Matrix 

The residual, V which is the difference between the 
estimate and the observed is usually useful when applying 
least squares adjustment technique for the determination of 
geoid heights interpolation model parameters since it is 
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equal to the difference between the model geoid heights 
and the computed geoid heights from GNSS observations 
and spirit levelling of the points. So, it can be used as a 
check. The unknown/model parameter is computed as 
(Oluyori et al., 2019) 

X = (A T A)~ l A T L (18) 

Where, 

(A T A) 1 = Inverse of the normal matrix 


are used for the computation of the RMSE, as well as the 
accuracy of the model. The Root Mean Square Error, 
RMSE index for the computation of the gravimetric- 
geometric geoid model accuracy as given by Yilmaz and 
Kozlu (2018) is 

RMSE = ±1— ^ (<W RcsidU aj ) 2 (19) 

V n i= i 

Where, 


The step by step procedures for the computation of the 
polynomial interpolation model Coefficients (CL mn ) of the 

variables (x, y) are detailed in Eteje and Oduyebo (2018). 
1.12 Accuracy of the Model 

The accuracy of the gravimetric-geometric geoid 
models is obtained using the Root Mean Square Error, 
RMSE index. To compute the accuracy of the model, the 
geoid heights of the selected points from the geoid model 
are compared with their corresponding gravimetric geoid 
heights to obtain the geoid height residuals. The geoid 
height residuals and the total number of the selected points 


/5V =N -N j 

Ui v Residual 1 v Computed 1 v Model ' 

N computed = Gravimetric Geoid Height 

Model = Gravimetric-Geometric Geoid Height 
n = Number of Points 

II. METHODOLOGY 

The methodology was divided into data acquisition, 
data processing and results presentation and analysis. 
Figure 3 shows the flow chart of the adopted methodology. 



Fig 3: Flow Chart of the Adopted Methodology 


2.1 Data Acquisition 

A total of 49 points were used in the study. The points 
included two primary control stations (XSU 92 and XSU 
100 were respectively located in Edo College and School 
of Nursing premises). The points were selected along the 
major roads of the City (See Figure 4). GNSS observation 


was carried out using CHC 900 dual-frequency GNSS 
receivers to obtain the coordinates and ellipsoidal heights 
of the points. The observations were carried out relative to 
control station XSU 92 using the static method of GNSS 
observation (See Figures 5 and 6). 
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Fig 4: Selected GPS and Gravity Points 


Fig. 5: Base Receiver at Control 
Station XSU92 


Fig. 6: Rover Receiver at One of the 
Points (RR01) at Ring Road 




The selected points were observed with a gravimeter 
(SCINTREX CG-5 Autograv Gravimeter) to obtain their 
absolute gravity values. The observation was carried out 
by an expert, a Geophysicist from Mountain Top 
University, Ibafo, Ogun State. The gravity observation of 


the points was carried out in seven different loops relative 
to a point whose absolute gravity value was known which 
was located within the Benin City Airport premises (See 
Figures 7 and 8). 



Fig. 7: Gravimeter Set over Reference Station at Benin City Airport Fig. 8: Gravimeter at Control Station XSU92 


2.2 Data Processing 

The GNSS observations were respectively downloaded 
and processed with HcLoader and Compass Post¬ 
processing software to obtain the positions and the 
ellipsoidal heights of the points. The geographic and 
rectangular coordinates, as well as the ellipsoidal heights 
of the points, were processed in Minna datum. The gravity 
observations of the points were processed by the expert 
who carried out the observations to obtain their absolute 
gravity values. All the necessary corrections such as drift 
correction, etc were applied during the processing. The 


theoretical gravity values of the points were computed on 
the local (Minna) datum ellipsoid (Clarke 1880 ellipsoid) 
using the latitude coordinates of the points, as well as 
equation (3). The gravity anomalies of the points were 
computed by finding the differences between the absolute 
gravity values of the points and their corresponding 
theoretical gravity values, as well as using equation (4). 
The computation of the free air correction requires the 
application of the orthometric heights of the points. And 
these were obtained by interpolation using the orthometric 
heights and the absolute gravity values of the two primary 
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control stations (XSU 100 and XSU 96). The orthometric 
heights of the points were interpolated as there was no 
Digital Terrain Model (DTM) of the study area. The free 
air correction was applied to the computed gravity 
anomalies of the points using equation (5). The free air and 
the Bouguer gravity anomalies of the points were 
computed but the free air gravity anomalies were used in 
the study. This is because the geoid heights of the two 
primary control stations obtained from their known 
orthometric and ellipsoidal heights approximated the geoid 
heights of the stations computed using the free-air gravity 
anomalies, as well as equation (6). The gravimetric geoid 
heights of the points were computed with the geographic 
coordinates, free air gravity anomalies and the theoretical 
gravity values of the points using equation (6). The 
computation of the gravimetric geoid heights of the points 
required the application of the spherical radius and it was 
computed using equation (7). The computed gravimetric 
geoid heights of the points using equation (6) were co- 
geoid heights. To obtain precise gravimetric geoid heights 
of the points, the combined topographic effect has to be 
computed and applied to the co-geoid heights. The 
combined topographic effect was computed using equation 
(8). To obtain the best gravimetric-geometric geoid model 
of Benin City, three geometric geoid surfaces (plane 
surface, equation (9), second degree surface, equation (10) 
and third degree surface, equation (11)) were fitted to the 
gravimetric geoid heights of the points. The fitting of the 
geometric geoid surfaces to gravimetric geoid heights of 
the points requires the computation of the model 
parameters using the least squares adjustment technique. 
The models' parameters were computed using equation 
(18). Microsoft Excel programs were developed with the 
computed models' parameters. The RMSE, as well as the 
accuracy of the models, were respectively computed with 
equations (19). 


III. RESULTS PRESENTATION AND 
ANALYSIS 

3.1 Analysis of the GNSS Observation Results 

The DGPS observation was carried out to obtain the 
coordinates and the ellipsoidal heights of the selected 
points. The DGPS observations were processed using 
Compass post-processing software. From the processing of 
the DGPS observations results, it was seen that the 
processed observations passed both the Network 
Adjustment Test and the % -Square (Chi-square) Test. 
This implies that the normal matrix generated was a 
regular one and inverted accordingly for the calculation of 
residuals. 

3.2 Analysis of the Gravimetric-Geometric Geoid 
Models 

Table 1 presents the gravimetric geoid heights; the 
three gravimetric-geometric geoid models (plane, second 
degree and third degree surfaces/models) heights and their 
respective RMSE while Figure 10 shows the plot of the 
RMSE of the three gravimetric-geometric geoid models. 
This was done to determine which of the three 
gravimetric-geometric geoid models is most suitable for 
application in the study area. The smaller the RMSE, the 
better the gravimetric-geometric geoid model. It is seen 
from Table 1 and Figure 10 that the RMSE of the plane 
surface, second degree surface, and the third degree 
surface models are respectively 0.9839m, 0.7126m and 
0.6746m which shows that the third degree gravimetric- 
geometric geoid model is most suitable for application in 
the study area. It can also be seen from Table 1 that the 
minimum and the maximum geoid heights from the third 
degree surface gravimetric-geometric geoid model are 
respectively 1.002m and 3.760m. This shows that geoid 
heights can be interpolated with the chosen (third degree 
surface gravimetric-geometric) geoid model within the 
range of 1.002m to 3.760m in Benin City, Nigeria. 


Table 1: Gravimetric Geoid Heights, the Three Gravimetric-geometric Geoid Models heights and their Absolute Differences 

between the Gravimetric Geoid heights 


POINT 

GRAVIMETRIC 

GEOID 

HEIGHT (m) (A) 

MODEL GEOID HEIGHT 

DIFF. BAY 

(A) & (B) 
(m) 

SQUARED 

DIFF. B/W 

(A) & (C) 
(m) 

SQUARED 

DIFF. B/W 

(A) & (D) 
(m) 

SQUARED 

PLANE 

SURFACE 
(m) (B) 

SECOND 

DEGREE 
(m) (C) 

THIRD 

DEGREE 

(m) (D) 

XSU92 

2.086 

2.386 

0.749 

1.072 

0.0900 

1.7876 

1.0282 

RR01 

2.420 

2.069 

1.456 

2.257 

0.1232 

0.9293 

0.0266 

SR01 

1.588 

1.989 

1.607 

2.275 

0.1608 

0.0004 

0.4720 

SR02 

1.978 

1.927 

1.791 

2.209 

0.0026 

0.0350 

0.0534 

SR04 

2.520 

1.799 

2.504 

3.260 

0.5198 

0.0003 

0.5476 

SR05 

2.802 

1.614 

3.266 

3.445 

1.4113 

0.2153 

0.4134 

SR06 

3.266 

1.560 

3.442 

2.961 

2.9104 

0.0310 

0.0930 

XSU100 

2.098 

1.897 

1.872 

2.214 

0.0404 

0.0511 

0.0135 
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AR01 

0.685 

1.922 

1.239 

1.240 

1.5302 

0.3069 

0.3080 

AR02 

0.720 

1.698 

1.343 

1.082 

0.9565 

0.3881 

0.1310 

AR03 

1.436 

1.425 

2.452 

2.233 

0.0001 

1.0323 

0.6352 

AR04 

1.439 

1.342 

2.827 

2.658 

0.0094 

1.9265 

1.4860 

UU01 

4.658 

2.307 

3.880 

3.667 

5.5272 

0.6053 

0.9821 

UU02 

3.498 

2.259 

2.527 

2.905 

1.5351 

0.9428 

0.3516 

UU03 

1.981 

2.197 

1.680 

1.761 

0.0467 

0.0906 

0.0484 

UU04 

1.276 

2.142 

1.538 

1.407 

0.7500 

0.0686 

0.0172 

UU05 

1.346 

2.208 

1.429 

1.330 

0.7430 

0.0069 

0.0003 

UU06 

1.489 

2.228 

0.940 

1.175 

0.5461 

0.3014 

0.0986 

UU07 

1.329 

2.122 

1.253 

1.197 

0.6288 

0.0058 

0.0174 

UU08 

1.263 

2.039 

1.469 

1.374 

0.6022 

0.0424 

0.0123 

AD01 

2.986 

2.777 

2.171 

2.299 

0.0437 

0.6642 

0.4720 

AD02 

4.019 

2.956 

2.839 

3.041 

1.1300 

1.3924 

0.9565 

AD03 

4.420 

3.059 

3.103 

3.321 

1.8523 

1.7345 

1.2078 

AK01 

1.376 

2.217 

1.211 

1.211 

0.7073 

0.0272 

0.0272 

AK02 

1.473 

2.360 

1.242 

1.002 

0.7868 

0.0534 

0.2218 

AK03 

2.252 

2.583 

1.503 

1.463 

0.1096 

0.5610 

0.6225 

AK04 

3.101 

2.816 

2.240 

2.485 

0.0812 

0.7413 

0.3795 

AK05 

3.954 

3.026 

3.615 

3.283 

0.8612 

0.1149 

0.4502 

MR01 

1.300 

2.191 

1.415 

1.372 

0.7939 

0.0132 

0.0052 

MR02 

1.488 

2.386 

1.560 

1.220 

0.8064 

0.0052 

0.0718 

MR03 

1.614 

2.561 

2.115 

1.786 

0.8968 

0.2510 

0.0296 

MR04 

4.037 

2.919 

3.119 

3.146 

1.2499 

0.8427 

0.7939 

MR05 

4.313 

3.010 

3.069 

3.302 

1.6978 

1.5475 

1.0221 

SK01 

1.500 

2.163 

1.193 

1.043 

0.4396 

0.0942 

0.2088 

SK02 

2.035 

2.219 

1.542 

1.233 

0.0339 

0.2430 

0.6432 

SK03 

2.379 

2.275 

2.161 

1.931 

0.0108 

0.0475 

0.2007 

EKS 

0.665 

2.001 

1.204 

1.325 

1.7849 

0.2905 

0.4356 

SLKO 

0.668 

1.978 

1.222 

1.051 

1.7161 

0.3069 

0.1467 

SLK01 

0.781 

1.915 

1.345 

1.067 

1.2860 

0.3181 

0.0818 

SLK02 

1.326 

1.846 

1.942 

1.716 

0.2704 

0.3795 

0.1521 

SLK03 

1.736 

1.745 

2.537 

2.506 

0.0001 

0.6416 

0.5929 

SLK04 

2.688 

1.702 

3.314 

3.431 

0.9722 

0.3919 

0.5520 

SLK05 

3.357 

1.699 

3.555 

3.760 

2.7490 

0.0392 

0.1624 

EK01 

0.078 

1.757 

0.776 

1.043 

2.8190 

0.4872 

0.9312 

EK02 

0.983 

1.439 

1.643 

2.010 

0.2079 

0.4356 

1.0547 

EK03 

1.729 

1.124 

3.160 

3.156 

0.3660 

2.0478 

2.0363 

EK04 

2.001 

0.969 

3.389 

3.344 

1.0650 

1.9265 

1.8036 

EK05 

2.516 

0.701 

2.480 

2.247 

3.2942 

0.0013 

0.0724 

AIRPORT 

0.578 

1.705 

1.298 

1.058 

1.2701 

0.5184 

0.2304 

RMSE = 

0.9839 

0.7126 

0.6746 
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1.2 

1 

OS 

0.6 

0.4 

0.2 

0 

Third Second Plane 

Degree Degree Surface 

Fig 10: Plot ofRMSE of the three Gravimetric-Geometric 
Geo id Models 

Figures 11a to lid respectively show the surface plots of 
the gravimetric, third degree surface gravimetric - 
geometric, second degree surface gravimetric-geometric 
and the plane surface gravimetric-geometric geoid heights 
of the selected points. This was done to present 
graphically the shape of the gravimetric geoid heights and 
those of the three gravimetric-geometric geoid models to 
determine which of the three gravimetric-geometric geoid 
models shape is most identical with the gravimetric geoid 
heights shape by comparing the three gravimetric- 
geometric geoid models shapes with that of the 

gravimetric heights. Comparing the three gravimetric- 
geometric geoid models shapes with that of the 

gravimetric geoid heights also enables the determination 
of which of the three gravimetric-geometric geoid models 
is most suitable for application in the study area. It can be 
seen from Figures 11a to lid that the shape of the third 
degree surface gravimetric-geometric geoid model is most 
identical with that of the gravimetric geoid heights. This 
also implies that the third degree surface gravimetric- 
geometric geoid model is most suitable for application in 
the study area. 

i 




Fig 11a: Surface Plot of Gravimetric Geoid Heights 



Fig lib: Surface Plot of Third Degree Gravimetric- 
Geometric Model Geoid Heights 



Fig 11c: Surface Plot of Second Degree Gravimetric- 
Geometric Model Geoid Heights 



Fig lid: Plot of Plane Surface Gravimetric-Geometric 
Model Geoid Heights 

IV. CONCLUSIONS AND 
RECOMMENDATIONS 

1. The study has established the local geoid model of 
Benin City, Nigeria from the combination of 
gravimetric and geometric methods of geoid 
determination. 

2. The determined geoid model was termed local geoid 
model as the absolute gravity values used were 
obtained from gravimetric observations carried out 
physically in the study area and as their respective 
theoretical gravity values were computed on the 
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local ellipsoid adopted for geodetic computation in 
the study area. 

3. The gravimetric-geometric geoid models were 
obtained by fitting three geometric geoid surfaces to 
the computed gravimetric geoid heights using the 
least squares technique. 

4. The accuracy of the three gravimetric-geometric 
geoid models was obtained by computing their 
respective RMSE. 

5. The third degree surface gravimetric-geometric 
geoid model is recommended for application in 
Benin City, Nigeria as it has the highest accuracy 
among the three gravimetric-geometric geoid 
models. 

6. A Microsoft Excel program was developed for full 
implementation of the determined geoid model in 
the study area. 

7. It is recommended that whenever ellipsoidal heights 
obtained from GNSS observations within Benin 
City, Nigeria are to be converted to practical, as well 
as orthometric heights, the determined geoid model 
should be applied. 

8. It is also recommended that the use of assumed, as 
well as handheld GPS receiver heights for 
engineering works should be totally abolished as this 
study has established the local geoid model of Benin 
City. 
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